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Abstract 

We study the coupled translational, electronic, and field dynamics of the combined 
system "a two-level atom + a single-mode quantized field + a standing-wave ideal 
cavity". In the semiclassical approximation with a point-like atom, interacting with 
the classical field, the dynamics is described by the Heisenberg equations for the 
atomic and field expectation values which are known to produce semiclassical chaos 
under appropriate conditions. We derive Hamilton - Schrodinger equations for prob- 
ability amplitudes and averaged position and momentum of a point-like atom in- 
teracting with the quantized field in a standing-wave cavity. They constitute, in 
general, an infinite-dimensional set of equations with an infinite number of integrals 
of motion which may be reduced to a dynamical system with four degrees of free- 
dom if the quantized field is supposed to be initially prepared in a Fock state. This 
system is found to produce semiquantum chaos with positive values of the maximal 
Lyapunov exponent. At exact resonance, the semiquantum dynamics is regular. At 
large values of detuning \6\ » 1, the Rabi atomic oscillations are usually shallow, 
and the dynamics is found to be almost regular. The Doppler - Rabi resonance, 
deep Rabi oscillations that may occur at any large value of \6\ to be equal to \apo\, 
is found numerically and described analytically (with a to be the normalized recoil 
frequency and pq the initial atomic momentum). Two gedanken experiments are 
proposed to detect manifestations of semiquantum chaos in real experiments. It is 
shown that in the chaotic regime values of the population inversion ^^outj measured 
with atoms after transversing a cavity, are so sensitive to small changes in the initial 
inversion zjn that the probability of detecting any value of Zout in the admissible 
interval [—1, 1] becomes almost unity in a short time. Chaotic wandering of a two- 
level atom in a quantized Fock field is shown to be fractal. Fractal-like structures, 
typical with chaotic scattering, are numerically found in the dependence of the time 
of exit of atoms from the cavity on their initial momenta. 
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1 Introduction 



The correspondence between quantum and classical worlds has been a subject 
of much interest from the early days of quantum mechanics. The emergence 
of classical dynamical chaos from quantum mechanics is one of the most de- 
bated problem in this field [1]. Isolated and bounded quantum systems do not 
show sensitive dependence on initial conditions in the same way as classical 
systems due to discreteness of quantum phase space and lacking of notion of 
trajectories in quantum mechanics. The Schrodinger equation for an isolated 
quantum system demonstrates only (quasi) periodic solutions even if the classi- 
cal counterpart of the quantum system under consideration would be chaotic. 
However, it is valid only if the quantum system is assumed to be absolutely 
isolated from the surroundings. 

In this paper, we study the temporal evolution of the strongly coupled atom- 
field system consisting of a single two-level atom interacting with a single 
mode of the quantized field in a standing-wave ideal cavity without any leak- 
age of photons. If the atom is assumed to be at rest the respective Jaynes - 
Cummings Hamiltonian is known to be integrable in the rotating-wave ap- 
proximation both under fully quantum and semiclassical descriptions [2]. If 
the atom moves with a constant velocity (the Raman - Nath approxima- 
tion) in the direction along which the coefficient of the atom-field coupling 
is changed (say, periodically) the semiclassical evolution of the expectation 
values of the atomic and field operators has been shown to be chaotic [3,4,5] 
with positive values of the maximal Lyapunov exponent in respective ranges 
of control parameters. The same time-periodic Hamiltonian has been shown to 
generate quasiperiodic solutions of the respective time-dependent Schrodinger 
equation for the probability amplitudes [6]. The Jaynes - Cummings system 
possesses two degrees of freedom, the electronic (internal) atomic one and 
the field one. In fact, when emitting and absorbing photons, an atom not 
only changes its internal state but its velocity is changed as well due to the 
photon recoil. It is a pronouncing effect with cold atoms in a laser field (see 
for a review, for example, [7]). Taking into account the translational (exter- 
nal) atomic degree of freedom, we get the autonomous Hamiltonian (1) with 
three degrees of freedom. The respective semiclassical equations of motion for 
the expectation values of the atomic position and momentum operators, the 
atomic population operator, and the combined atom-field operators have been 
shown to be chaotic [8,9] with positive values of the maximal Lyapunov expo- 
nent. The semiclassical atom-field dynamics has been shown to demonstrate 
many interesting features, including the interaction of nonlinear resonances 
[10], atomic fractals [5,11], Levy flights and anomalous atomic diffusion [12], 
and the Doppler - Rabi resonance [13]. 
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In this paper, we go further in quantizing the atom-photon interaction in 
a standing-wave cavity. The field and internal atomic degrees of freedoms 
are treated as fully quantum ones obeying the Schrodinger equation. They 
are coupled to the external atomic degree of freedom obeying the Hamilton 
equations. Such a quantum-classical hybrid may be considered as a reality- 
based model of interaction between a quantum system and the surroundings. 
From the standpoint of dynamical system theory, the hybrid is described 
by an infinite-dimensional set of nonlinear ODE's with an infinite number 
of integrals of motion with clear physical meanings (see Eqs. (18) - (21)). 
Any state of the quantized field may be represented as a superposition of a 
number of the so-called Fock states where n is the number of photons 
in the respective state. Any state of a two-level atom is a superposition of 
its ground and excited states. The Hilbert space of the quantized atom-field 
subsystem (which should be treated as a whole unity) is an infinite direct 
sum of two-dimensional subspaces in each of which the so-called number of 
excitations should be conserved in the process of evolution. The atom-field 
quantized subsystem evolves in such a way that transitions, belonging to the 
subspaces with different values of the number of excitations, are forbidden. 
The main aim of the paper is to investigate the effects of quantization on those 
properties of the system that produce dynamical chaos. 

This paper is organized as follows. In Sec. 2 we introduce the system under 
consideration and the respective Hamiltonian. In Sec. 3 we derive the Heisen- 
berg equations for the expectation values of the atomic and field operators and 
review briefly the properties of semiclassical chaos. Our main results are given 
in Sec. 4 where we derive the Hamilton - Schrodinger equations of motion 
for a two-level atom in a quantized field, study the quantum Doppler - Rabi 
resonance, the properties of semiquantum chaos, and atomic fractals in the 
Fock quantized field. 



2 Two-level atom with recoil in a standing-wave cavity 



We consider a single two-level atom with the frequency oOa of an electric dipole 
transition and mass rUa moving in an ideal cavity which sustains a single 
standing-wave mode along the axis x with the frequency ujf and the wave 
vector kf. In the strong coupling limit, where the coefficient of the coupling 
fig is much greater than all the relaxation rates, the atom-field dynamics may 
be treated as Hamiltonian with the respective operator 

H — 1 — — ^ cTz + tvjjfa)a — HHq (a' a- + da+) coskfX, (1) 
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whose summands describe the kinetic and internal energies of the atom, the 
field energy, and the energy of the atom-field interaction, respectively. The 
momentum P, position X, atomic a, and field a, operators satisfy the 
standard commutations relations: 

[X, P] = ih, [a±, 0-^] = t2ct±, [a+, = a^, [a, d^ = 1. (2) 

Operators belonging to different degrees of freedom commute with each other 
at the same time moment. 

In the process of emitting and absorbing photons an atom not only changes its 
internal electronic state but its external translational state is changed as well 
due to the photon recoil effect. An interplay between the electronic, transla- 
tional, and field degrees of freedom of the strongly coupled atom-field system 
may be described as in the Heisenberg as in the Schrodinger pictures. 



3 Heisenberg equations for the atomic and field expectation values 
and semiclassical dynamics 



It is convenient to write down the Heisenberg equations for the following 
operators: 

- P ^ d'^a^ + «a+ aJa_ — da^ 
x^kfX, p = —, u = ^= , v = i j= , a^, (3) 

y N y N 



where = dVx + (a^ -|- /)/2 is a constant operator of the total number of 
excitations and / the identity operator. The derivative of an arbitrary operator 
A with respect to the normalized time r = ^ot 

ihA = [A, H] (4) 

results in the following Heisenberg equations for the operators (3): 

X — ap, 

p — — VA" usinx, 

u — Si), (5) 
V — — 5u + 2\f^ az cos X, 
Gz— — 2\/7^'f)cosx, 
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where the control parameters 



— ^^n 

are the normalized recoil frequency and the detuning between the field and 
the atomic frequencies, respectively. The parameter a = 2uir/^o characterizes 
the average change in translational energy tvujR in the process of emission and 
absorption of a photon. The set (5) is not closed. In order to describe fully 
quantized dynamics one should write down the equations of motion for all the 
products of the operators (3) and their functions which, in turn, would gener- 
ate another operator products and respective equations of motion and so on. 
This process, in general, generates an infinite hierarchy of operator equations. 
The simplest way to resolve this problem is to take quantum expectation val- 
ues over an initial quantum state and to factorize all the operator products in 
Eqs. (5). In this way one gets the closed dynamical system 

X — ap, 

p = —\fN 14 sin X, 

ii = Sv, (7) 
V = —Su + 2\fN z cos a;, 
z = —2\fN V cos a;, 

for the classical variables, namely, the atomic position x = (x) and momen- 
tum p = (p), the atom-field variables u = (ii,) and v = {v), and the atomic 
population inversion z = (cr^). A conserved number of excitations in the sys- 
tem = (N) = n -|- (2; -|- 1)/2 is the additional control parameter. The set (7) 
has two integrals of motion 

W ^^p'^ -uVNcosx - -z, i?^ = + + (8) 

where W is the conserved total energy and reflects the conservation of 
the length of the Bloch vector in the limit of the large number of photons, 
n — {a^d) ^ 1. The set (7) (with slightly another normalization) has been 
derived in [8,9]. 

In deriving Eqs. (7) from an infinite hierarchy of Heisenberg operator equa- 
tions, we treat an atom as a point particle that may be justified if its mo- 
mentum is much greater then the photon momentum Hkf, i. e. if \p\ ^ 1. 
Factorization of all the operator products in the respective operator equations 
means that we do not take into account either quantum nature of the field 
or quantum correlations between all the atomic and field degrees of freedom. 
It is justified if n 3> 1. Moreover, the procedure of reducing operator equa- 
tions to semiclassical ones is not unique because the form of the resulting 
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semiclassical equations depends on the factorization procedure. Semiclassical 
equations of motion, different from (7), have been derived in [5,12] with the 
same Hamiltonian (1). All of them, of course, has the same form in the limit 
n > 1. 

It has been found in a series of papers [8,9,5,10,12] that the semiclassical 
equations of motion of the strongly coupled atom-field system in a standing 
wave cavity produce different types of motion. At exact resonance, 6 = 0, 
the motion is regular since the set (7) gains an additional conserved quantity 
u{t) = u{0) = uq which reflects the conservation of the atom-field interaction 
energy at 5 = 0. Depending on the values of the initial atomic momentum 
Po, an atom either oscillates periodically in a potential well of the stand- 
ing wave or flies over its potential hills. In resonance, the optical potential 
U = —Uq^/N cosx — Sz/2 coincides with the standing wave structure. The 
oscillations of the internal atomic energy, the so-called Rabi oscillations z{t), 
and the oscillations of the atom-field variables, u{t) and v{t), are regular as 
well. The center-of-mass motion of the atom does not depend on the Rabi os- 
cillations, but its frequency depends on the initial interaction energy Uq since 
it determines the depth of the optical potential wells. In contrary, the Rabi 
oscillations depend on the translational motion since the strength of the atom- 
field coupling depends on the position of an atom in a cavity. The respective 
exact solutions of Eqs. (7) with S = one can found in [9]. 

Out of resonance, S 0, the set (7) with two integrals (8) is an autonomous 
Hamiltonian system whose motion takes place on a three-dimensional hyper- 
surface. It has, generally speaking, a positive Lyapunov exponent A which has 
been computed in the paper [8] as a function of the control parameters a and S 
and of the initial atomic momentum po- It follows from the first two equations 
in (7) that the atomic center-of-mass motion is described by the equation for 
a physical pendulum with a frequency modulation 

X + aVN u{t) sin x = 0, (9) 

where m is a function of time and all the other variables. Following to [14,10] 
one can show that Eq. (9) may produce a stochastic layer in a neighbourhood 
of the unperturped separatrix. Analogously to what has been done in [10], 
the normalized (to a separatrix value) width of the stochastic layer may be 
estimated as follows 



A ~ 87r(l]/a;)^exp ( — , (10) 



7in\ 



where = \/S^ + AN is the normalized Rabi frequency and the frequency 
uj = \j2aN^/^\6\/il characterizes small-amplitude oscillations. It should be 
noted that A gives the lower bound for the layer width (see [10]). 
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Center-of-mass motion of a two-level atom in an ideal standing-wave cavity 
has been found in [12] to be anomalous. A typical chaotic atomic trajectory 
consists of intervals of regular motion with an almost constant velocity in 
each interval (Levy flights) interrupted by erratic walks. Such an intermit- 
tency is typical of Hamiltonian systems with nonhomogeneous phase space 
with a fractal-like structure consisting of KAM tori, cantori, chains of islands, 
stochastic sea, etc. [15]. From this point of view, the Levy flights may be 
understood as those trajectories that "stick" to islands boundaries for a long 
time. The representative point on a typical chaotic trajectory sooner or later 
approaches, as closely as desirable, an island boundary that separates regular 
and chaotic motions. Nearby such a boundary, the maximal Lyapunov expo- 
nent A goes to zero, and cantories block the trajectory escape to the stochastic 
sea. As a result, the atomic motion is almost regular for a time that may be 
very long. From the physical point of view, the intermittency of the center- 
of-mass motion is due to the intermittent oscillations of the effective optical 
potential C/eff = —y/N ucosx — 5z/2 which governs the translational motion 
of the atom [13]. 

The Levy flights impact the statistical properties of the atomic motion result- 
ing in the anomalous diffusion. It was numerically found in [12] that the second 
moment of the position of atom in a cavity evolves in time as ~ r'*, where 
the transport exponent (i may vary from the value /i ~ 1 . corresponding to the 
normal diffusion, to the value /i ~ 2 corresponding to the superdiffusion. The 
Poincare theorem states that every trajectory of a closed conservative dynami- 
cal system, except for trajectories of the set of zero measure, returns arbitrarily 
close to its origin infinitely many times. The recurrence time distribution in 
the system with perfect mixing is known to be Poissonian, -P(t) = h~^e~^'^ , 
where h is the Kolmogorov-Sinai entropy. The motion with intermittency and 
Levy flights leads to the power law, -P(t) ~ r"''' at r — > cxo. The exponents ii 
and 7 are related to each other, and their values depends on the values of the 
control parameters a, 5, and N because by changing them one changes the 
topology of the phase space. 

The atom-photon interaction in a cavity may be considered as the chaotic 
scattering problem [5,11] where a two-level atom is scattered by the standing- 
wave light or by an optical potential. In difference from the usual scattering of 
atoms by light [7], the effective optical potential Ucs of the strongly coupled 
atom-field system depends not only on the field but on the atomic variables 
as well. Let two atomic detectors to be placed at the cavity mirrors and let 
them detect the time T of atomic exit from the cavity. Let the identically 
prepared atoms with given initial momenta po be placed one by one in the 
middle of the cavity. The dependence T{pq) has been found in [5,11] to have 
a beautiful selfsimilar structure with the Hausdorff dimension to be equal to 
d ~ 1.84. Tiny interplay between all the degrees of freedom is responsible for 
trapping atoms with T — > oo even in a very short microcavity. Simulation in 
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the cited papers has been performed with a cavity whose length is equal to two 
standing-wave lengths. Two kinds of atomic fractals have been found in [11], a 
countable fractal (a set of po generating separatrix-like atomic trajectories) and 
a seemingly uncountable fractal with a set of po generating infinite walkings 
of atoms inside the cavity. 



4 Hamilton-Schrodinger equations and semiquantum dynamics 

4.1 Derivation of equations of motion 

In this section we again consider a two-level atom as a point particle but 
electronic and field degrees of freedom are treated as fully quantum ones in the 
Schrodinger picture. It enables us to study the role that field quantum statistic 
and atomic superposition play in the full atom-field dynamics including chaos. 
We start with the Hamiltonian H (1). 

The Hamilton equations for the classical translational degree of freedom is 
easily found 

d{X)_^d{H)_ d{P)_^_d{H)_ 

dt d{p)' dt d{xy ^ ' 

where (...) denotes an expectation value of the corresponding operator over 
a quantum state l^') of the electronic-field Hamiltonian. Using the same nor- 
malizations and notations as in the preceding section, we get 

X = ap, p = —{u)smx, (12) 

where (u) = (^(r)|Mo|^(r)). 

Let us expand a state vector of the electronic-field subsystem over the basic 
energetic atomic states |2) and |1) and the Fock field states \n) 

= E(«n(T)|2, n) + 6„(t)|1, n)), (13) 

n=0 ^ ^ 

where a„ and 6„ are the probability amplitudes to find the atom in its excited 
or ground state with n photons in the mode, respectively. Substitution of the 
vector (13) in the time-dependent Schrodinger equation 

i^|^) = ^|^) (14) 
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gives the infinite-dimensional set of the coupled ODE's 

d„ = -i [Aaan - Vn+ 1 fen+i cosx) , 

; ^ { (15) 

bn+i= i Ab6* - + 1 a* COS x j , n = 0, 1, 2, . . . , 

where Aa = ncuf + u!a/'2 and A;, = {n + l)u!f — cua/'^- Introducing the following 
combinations of the probability amplitudes: 

M„ = 2Re (a„6;^_^i), v„ = -2 Im (an6*+i) , = |a„p - (16) 

we get the quantum Bloch-like equations 

Vn — —6Un + 2\/ n+ IZn COS X, (17) 
Zn — —2y/n+lVnCOSX, u — 0, 1, 2, ... , 

where the detuning S is the same as in (6). After computing {u), we obtain 
our basic Hamilton-Schrodinger equations 

X — ap, 

•p = — V JT- + 1 tin sin X, 

n=0 

Un = SVn, ^^^^ 
Vn = —SUn + 2\Jn + 1 Zn COSX, 



Zn = —2^/n + l Vn COS X , n = 0, 1, 2, . . . . 

This infinite set of nonlinear ODE's possesses an infinite number of the inte- 
grals of motion, the total energy integral 

a °° S °° 

= -p - ^ Vn + l Mucosa; -- (19) 

^ n=0 ^ n=0 

the Bloch-like integrals for each n 

Rl-< + < + 4, (20) 

and the global integral reflecting conservation of the total probability 

oo 

E = 1. (21) 
n=0 

The quantities that can be measured in real experiments are the atomic posi- 
tion X, momentum p and the atomic population inversion 

4r) = f:zn{r). (22) 

n=0 
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The semiquantum equations (18) — (21) should be compared with the semi- 
classical equations (7) and (8) because they describe the same physical situ- 
ation but on the different ground. They are exactly identical in the case of 
the initial Fock state of the cavity field with n quanta and the atom to be 
prepared initially in one of its energetic states. It is easy to show that an in- 
finite number of equations (18) reduces in this case to five equations (7) with 
N = n + {z + l)/2 = n+l if |^(0)) = \2, n) and N = n if |^(0)) = |1, n). If 
the atom is initially prepared in a general superposition state and the field is 
in the Fock state |n) 



l^'(O)) = a„(0)|2, n) + 6„(0)|1, n), a„(0) + 6„(0) = 1, 
we get from Eqs. (18) 



(23) 



X — ap, 

p = — (y^Un-i + Vn + 1 Mn) sinx, 

iln-l = 5Vn-l, Un = SVn, 

Vn-1 — —SUn-1 + 2\/^Zn-l COSX, = + 2\/n + 1 Zn COSX, 

Zn-1 — —2\/Evn-iCOSX, Zn = ~2^ n^\Vn COS X 



(24) 



with the respective integrals 



a 2 

W — —p — {^/nUr, 

= + + MO) 



„_i Vn + 1m„) cos a; - -{zn-\ + ^n), 



2 I 2 I 2 



Rn-1 + Rn — 1) 

and initial conditions 



an(0) 



(25) 



a;(0) = xo, p{0) = po, Zn-i{0) = - 6„(0) , 2;„(0) = a„(0) 

«„_i(0) = «„(0) = ^;„_i(0) = ^;„(0) = 0. (26) 



4-2 Doppler - Rabi resonance 



If the field frequency is far detuned from the frequency of the atomic 
working transition uja, i- e. if \5\ ^ 1, then the Rabi oscillations z{t) should 
be very shallow as oscillations of an oscillator under influence of a far-detuned 
time-dependent force. It is correct if the atom does not move. The Doppler 
effect with a moving atom causes an interesting effect of the Doppler - Rabi 
resonance if the condition \ap\ ^ \6\ is fulfilled. A standing wave is a sum of 
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two running waves moving in the opposite directions. In the reference frame 
of the moving atom, their frequencies are different due to the Doppler effect 

Va Va 

uji=u}f ujf, u}2 = ulf^ ujf, (27) 

where Va and c are the atomic and light velocities, respectively, and Va <^ c. 
If the initial atomic momentum po is sufficiently large, the Raman - Nath 
condition, p ~ po is valid. Let us define dimensionless detunings between the 
frequencies of the atomic transition and the running waves as follows: 

LOi-LOa r. r. ixl2- UJa r . fr,o\ 

di = — = d - apo, 62 = — = d + apQ. (28) 

It is evident that the atom comes in resonance with one of the running waves 
if |apo| — 1^1- If 1^1 3> 1, the interaction of the atom with the other running 
wave is negligibly small. Suppose that the atom is initially prepared in the 
ground state, i. e. Zn-i{^) — —1, -Zn(O) = 0, and z{Q) — —1, than the equations 
of motion (24) for the atom with a constant speed and with the other initial 
conditions written in (26) are reduced to the following simple set: 

Vn ^ -5Un + VnTT Zn, (29) 

Zn = -y/n+lVn- 

It should be noted that the amplitude of the running wave is half of the 
standing-wave amplitude, and the atom-field interaction does not depend on 
atomic position. The solution for the atomic population inversion with given 
initial conditions is easy to find 



z{t) = Znir) = - [ '^ Q^ j cosQ„T, (30) 



where n„ = y (5 — apo)^ + + 1 is the Rabi frequency. In particular, at 
exact Doppler - Rabi resonance, \apo\ = \6\, the inversion oscillates at the 
frequency {n+ 1)^/^ and its amplitude is maximal. 

The Doppler - Rabi resonance is illustrated in Fig. 1 where we plot the Rabi 
oscillation signals z{t) computed with the full set (24) and two initial states 
of the atom when it is prepared in the ground level with z{0) = —1 and in the 
superposition state with z{0) = {zn-i{0) = —1/2, and Zn{0) = 1/2). In our 
simulation a — 10~^, ri = 10 and xq — 0. When 5 — 32 andpo = 32000 and the 
resonance condition is fulfilled we really see the maximal Rabi oscillations with 
both the initial conditions. The signal in Fig. la with z{0) = —1 is described 
by the analytical solution (30). When z{0) = we have two oscillators with 
different frequencies (see (24)) and the resulting signal z{t) = 2;„_i(r) + Zn{T) 
in Fig. Ic demonstrates a beating with two harmonic components of the type 



11 




Fig. 1. Doppler - Rabi resonance at the condition apo = 5 with 5 = 32, po = 32000, 
a = 0.001, (a) z{0) = — 1 and (c) z(0) = 0. Shallow Rabi oscillations with the same 
value of po but out off the resonance are shown for comparison with (b) z{0) = — 1, 
(5 = 1 and (d) ^(0) = 0,6= 10. 

of (30) with Qn = {n + 1)^/^ and fin-i — n^^'^- For comparison, we plot in 
Fig. lb and d very shallow Rabi oscillations with the same value of the atomic 
momentum po = 32000 but out off the resonance, japol 7^ l*^!- The amplitude of 
the Rabi oscillations in Fig. lb with z{{)) = —1 and 5 = 1 does not exceed 2% 
of the amplitude of the resonant oscillations in Fig. la. The same is valid with 
another initial value of the atomic inversion population, 2(0) = (compare, 
please. Figs. Ic and d). 

In conclusion, deep Rabi oscillations are possible at as large values of the 
detuning \5\ as desirable if a two-level atom moves with the corresponding 
velocity. 



4-3 Chaos in a quantized Fock field 

Lyapunov exponents are known to be quantitative indicators of chaos in dy- 
namical systems. They characterize the behavior of close trajectories in phase 
space. If the quantized field is initially prepared in a Fock state \n) with n 
quanta in the mode, the Hamilton - Schrodinger equations of motion have the 
form (24) for an arbitrary internal atomic state. The 8-dimensional dynamical 
system (24) with 4 integrals of motion (25) has, as maximum, four nonzero 
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Lyapunov exponents Xi {i — 1, 2, 3, 4) 



A,; 



lim A, 

T— +0O 



[r), A,(r) 



lim — In . ^ ^, } , 

A,(0)-.OT Ai(0)' 



(31) 



where Aj(r) is the distance (in the Euclidean sense) in the i-th direction 
at the moment r between two trajectories that were close to each other at 
T = 0. In Hamiltonian systems, due to the phase space volume conservation, 
Ai + A2 + A3 + A4 = 0, and Ai = — A2, A3 = — A4. Computing the maximal 
Lyapunov exponent A, one measures an averaged rate of separation of initially 
close trajectories. 




Fig. 2. The maximal Lyapunov exponent A with zero initial atomic population inver- 
sion z(0) = and n = 10 photons in the Fock quantized field versus the atom-field 
detuning 6 and the logarithm of the dimensionless recoil frequency a. 

The system has three control parameters, the initial number of photons in 
the mode n, the normalized recoil frequency a and the detuning S, each of 
which may vary in a wide range of values. To diagnose chaos it is instructive 
to compute the so-called topographic A-maps [3,4,17] which show by color 
modulation values of A in the ranges of values of two control parameters with 
the other to be fixed. When computing the A-maps, we choose the following 



initial conditions: a;(0) = 0,p{0) = 50, m„_i(0) = m„(0) 



Jn-l 



(0) = v,,{0) = 0, 



Zn-i{0) = -1/2, Zn{0) = 1/2 and 2;„_i(0) = 0, 2;„(0) = 1. It means that the 
atom is prepared initially in the state with zero population inversion z{0) = 
or in the excited state, z{0) — 1. In Fig. 2 we show the A-map with z{0) = in 
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Fig. 3. The double logarithmic plot of A with an initially excited atom ^;(0) = 1 and 
n = 10 photons in the Fock quantized field versus a and n with 6 = 0.5. 

the ranges of the detuning |5| ^ 2 and the recoil frequency —4 ^ logio^ ^ ~1 
with the fixed value of the initial number of photons n — 10. The set (24) is 
integrable at exact atom-field resonance {S = 0) with A = 0. The a-S map 
confirms this conjecture. The values of a ~ 10~^ 10~^, that correspond to 
positive values of A, are reasonable with real atoms [9], and one may expect 
chaotic atomic motion in the respective ranges of a and 6. Another A-map 
with z{0) = 1 shows in Fig. 3 the value of A in dependence on a and n at 
S — 0.5 in double logarithmic scale. The magnitude of A grows, in average, 
with increasing the initial number of photons in the cavity mode. 

A feasible scheme for detecting manifestation of chaos with hot two-level Ry- 
dberg atoms moving in a high-Q microwave cavity has been proposed in [5]. 
The same idea could be realized with cold usual atoms in a high-Q micro- 
cavity. Consider a 2D-geometry of a gedanken experiment shown in Fig. 4a, 
where a monokinetic atomic beam propagates almost perpendicularly to the 
cavity axis x. In a reference frame moving with a constant velocity in the 
y-direction, there remains only the transverse atomic motion along the axis x. 
One measures atomic population inversion after passing the interaction zone. 
Before injecting atoms in the cavity, it is necessary to prepare all the atoms 
in the same electronic state, say, in the excited state, with the help of a vr- 
pulse of the laser radiation. It may be done with only a finite accuracy, say, 
equal to Azin for the initial population inversion Zi^. The values of the pop- 
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Fig. 4. Schematic diagrams showing scattering of atoms at the standing wave, (a) A 
monokinetic beam of atoms propagates transversally to the cavity axis x, and the 
atoms are detected outside the cavity, (b) Atoms are placed one by one inside the 
cavity with almost zero transversal velocity and are detected at the cavity mirrors. 

ulation inversion z^^i are measured with detectors at a fixed time moment. 
If we would work with the values of the control parameters corresponding 
to the regular atom-field dynamics, we would expect to have a regular curve 
^out^-2in- In the chaotic regime, the atomic inversion at the output can be pre- 
dicted (within a certain confidence interval Az) for a time not exceeding the 
so-called predictability horizon 



1 , Az 
m 



(32) 



which depends weakly on Az^^ and Az. Since the maximal confidence interval 
lies in the range \Az\ ^ 1 and A may reach the values of the order of 1.5 
(see A-maps), the predictability horizon in accordance with the formula (32) 
can be very short: with A =0.5 may be of the order of 10 in units of the 
reciprocal of the vacuum Rabi frequency VIq that correspond to tp ~ 10~^ s 
with the realistic value of VLq ~ 10*^ rad- s~^ [18]. 

In the regular regime, the inevitable errors in preparing Az-^-^ produce the 
output errors Az^^i of the same order. In the chaotic regime, the initial uncer- 
tainty increases exponentially resulting in a complete uncertainty of the de- 
tected population inversion in a reasonable time. It is demonstrated in Fig. 5, 
where we plot the dependence of the values of z{t) — Zout at r = 100 (Fig. 5a) 
and T = 200 (Fig. 5b) on the values of z{0) — Zi^ in the chaotic regime with 
6 = 0.4 and A ~ 0.05. It is evident from Fig. 5a, that at the detection time 
moment r = 100 an initial error Azi^ = 10^"^ leads to a complete uncertainty 
Azout — 2 in rather large vicinities of the initial inversion Zi^ ~ ±1, whereas 
the dependence -Zout^-^in is a regular one in the vicinity of — 0. With increas- 
ing the detection time (see Fig. 5b at r = 200), the probability of detecting 
any value of ^out in the interval [—1, 1] is almost unity in the whole range of 
Z[n. The formula (32) gives the value of the predictability horizon Tp ~ 200 
with A ~ 0.05. To feel the difference, it is desirable to carry out a control 
experiment at the exact resonance {S — 0) when the atomic motion is fully 
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Fig. 5. Dependence of the output values of the atomic population inversion Zout on 
its initial values Zin with S = 0.4 (a) at r = 100 and (b) at r = 200 with the inset 
showing this dependence at the exact atom-field resonance, 6 = 0. 

regular with any initial values. The dependence Zout^Zi^ with 5 = is demon- 
strated in the inset in Fig. 5b with all the other control parameters and initial 
values being the same. 



4-4 Atomic fractals in a quantized field 



In this section, we treat the atom-photon interaction in a high-Q cavity as 
a chaotic scattering problem [19,20]. Atoms from outside are injected into a 
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ID-cavity, interact for a while with the cavity field and then they are removed 
from the cavity by hook or by crook. Following to [5,11], let us consider the 
scheme of scattering of atoms by the standing wave shown in Fig. 4b. Atoms, 
one by one, are placed at the point x — with different initial values of the 
momentum po along the cavity axis. For simplicity, we suppose that they have 
no momentum in the other directions (ID-geometry). We compute the time 
the atoms need to reach one of the detectors placed at the cavity mirrors. 
The dependence of this exit time T on the initial atomic momentum Pq is 
studied under the other initial conditions and parameters being the same. To 
avoid complications that are not essential to the main theme of this section, 
we consider the cavity with only two standing-wave lengths. Before injecting 
into a cavity, atoms are suppose to be prepared in the superposition state with 
Un-i(O) = M„(0) = t'„-i(0) = vjO) = 0, 2„_i(0) = -l/2,z„(0) = 1/2, i. e. in 
the state with zero population inversion z{0) = Zn-i{0) + Zn{0) = 0. 

At exact resonance {5 = 0) with m„_i(t) = Unir) = 0, the optical potential 
U = (a/^m„_i + \/n + 1 Un) cosx — S{zn-i + Zn) /2 is equal to zero, and the 
analytical expression for the dependence in question can be easily found to 
be the following: T{6 = 0) = 37r/2apo if Po > and T{6 = 0) = n/2apo if 
Po < 0. Atoms simply fly through the cavity in one direction with their initial 
constant velocity and are registered by one of the detectors. Out of resonance 
(S 7^ 0), the atomic motion has been numerically found in the preceding 
section to be chaotic with positive values of the maximal Lyapunov exponent. 
Fig. 6 shows the function T(po) with the normalized detuning S = 0.4, the 
recoil frequency a — 10'^, and the average initial number of cavity photons 
n — 10. The exit-time function demonstrates an interniittency of smooth 
curves and complicated structures that cannot be resolved in principle, no 
matter how large the magnification factor. Fig. 6b shows magnification of the 
function for the small interval 64.1 ^ po ^ 64.6. Further magnification in 
the range 64.2743 ^ po ^ 64.2754 shown in Fig. 6c reveals a beautiful self- 
similar structure. Some structures in Fig. 6a that look like fractal are not, in 
fact, unresolvable and self-similar. Magnification of the structure in the range 
73.2 ^ Po ^ 73.8 (Fig. 7a), shown in Fig. 7b, demonstrates quite a smooth 
function without unresolvable substructures and with only two singular points 
on the borders of the respective momentum interval. Beating, visible in all 
the structures of the atomic fractal in Figs. 6, 7, should be attributed to 
the structure of the Hamilton-Schrodinger equations (24) which describe two 
atom-field oscillators with slightly different frequencies. 

The exit time T, corresponding to both smooth and unresolved po intervals, 
increases with increasing the magnification factor. It follows that there exist 
atoms never reaching the detectors in spite of the fact that they have no ob- 
vious energy restrictions to leave the cavity. Tiny interplay between chaotic 
external and internal dynamics prevents these atoms from leaving the cavity. 
The similar phenomenon in Hamiltonian systems is known as dynamical trap- 
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Fig. 6. Atomic fractal in a quantized Fock field with different resolutions {S = 0.4 
and z{0) = 0). 

ping [16]. Different kinds of atomic trajectories, which are computed with the 
system (24), are shown in Fig. 8. A trajectory with the number m transverses 
the central node of the standing-wave, before being detected, m times and is 
called m-th trajectory. There are also special separatrix-like m5'-trajectories 
following which atoms in infinite time reach the stationary points Xg — ±7rn 
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Fig. 7. Example of one of the non-fractal substructures in the atomic fractal shown 
in Fig. 6a with successive magnifications. 

(n = 0, 1, 2, . . . ), Ps = 0, transversing m times the central node. These points 
are the anti-nodes of the standing wave where the force acting on atoms is 
zero. A detuned atom can asymptotically reach one of the stationary points 
after transversing the central node m times. The trajectory wdth number 1, 
showing in Fig. 8, is close to a separatrix-like liS-trajectory. The smooth po in- 
tervals in the first-order structure in Fig. 6a correspond to atoms transversing 
once the central node and reaching the right detector. The unresolved singu- 
lar points in the first-order structure with T = oo at the border between the 
smooth and unresolved po intervals are generated by the IS'-trajectories. Anal- 
ogously, the smooth and unresolved po intervals in the second-order structure 
in Fig. 6b correspond to the 2-nd order and the other trajectories, respectively, 
with singular points between them corresponding to the 2S'-trajectories and 
so on. 



There are two different mechanisms of generation of infinite exit times, namely, 
dynamical trapping with infinite oscillations (m = oo) in a cavity and the 
separatrix-like motion (m 7^ 00). The set of all initial momenta generating 
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Fig. 8. Sample atomic trajectories transversing the central node in the scheme Fig. 4b 
m times (m = 1, 2, 3). The trajectory with number 1 is close to a separatrix-like 
l-S-trajectory. 

the separatrix-like trajectories is a countable fractal. Each point in the set can 
be specified as a vector in a Hilbert space with m integer nonzero components. 
One is able to prescribe to any unresolved interval of m-th order structure a set 
with m integers, where the first integer is a number of a second-order structure 
to which trajectory under consideration belongs in the first-order structure, 
the second integer is a number of a third-order structure in the second-order 
structure mentioned above, and so on. Such a number set is analogous to a di- 
rectory tree address: "<a subdirectory of the root directory>/<a subdirectory 
of the 2-nd level>/<a subdirectory of the 3-rd level>/...". 



5 Conclusion 

Atoms in cavities provide a microscopic nonlinear dynamical system with a 
reach variety of qualitatively different dynamics that may be explored to study 
the key problems of modern quantum physics and nonlinear science including 
the problem of the quantum-classical correspondence. The up-to-date experi- 
mental state of art has reached the stage where the quantum to classical transi- 
tion and the borderland between them can now be probed directly. Increasing 
the number of atoms or/and the average number of photons in a cavity mode, 
one can force the atom-field system to operating in quantum, semiquantum, 
and semiclassical regimes providing a link between micro-, meso-, and . In 
this paper, we derived the Hamilton - Schrodinger equations of motion that 
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describe the semiquantum Hamiltonian dynamics of a two-level atom with 
recoil strongly coupled to a single-mode standing-wave quantized field in an 
ideal cavity. It has been shown that dynamical chaos and fractals may arise in 
a wide range of the control parameters of the system. To estimate their mag- 
nitudes we use the parameters of the real experiments with single atoms in 
Fabry - Perot microcavities in the strong-coupling regime [18], for which the 
amplitude of the atom-field coupling strength may reach Qq = 271 ■ (10^^ 10®) 
Hz exceeding the decay rates of the cavity mode and the atomic dipole. With 
the above mentioned values of Qq and A;/ ~ 27r • 10^ m~^, one can estimate 
the normalized recoil frequency a to be in range 10~^-10~^ depending on the 
magnitude of the atomic mass. 
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